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Abstract 

It is shown that in extremely anisotropic layered superconductors the interac- 
tion of vortex lines with a parallel planar surface, which for straight lines along 
the c-axis decreases exponentially over the in-plane penetration depth A, be- 
comes a long-range dipole-dipole attraction when the vortex line is distorted 
randomly. This novel long-range fluctuation-induced attraction enhances the 
thermal fluctuations down to depths much larger than A and may lead to flux 

creep towards the surface. 

PACS numbers: 74.30. Gn, 74.60.Ge 
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Abrikosov vortex lines in layered superconductors have several unusual properties as 
compared to vortex lines in isotropic or weakly anisotropic superconductors. In particular, 
vortex lines oriented perpendicular to the superconducting layers may be considered as a 
stack of two-dimensional point vortices or pancakes |1]-§J. In the case of large anisotropy 
the pancakes interact via a magnetic pair-potential, which parallel to the layers decreases 
logarithmically and perpendicular to the layers decreases exponentially. Most importantly, 
the interaction of pancakes within the same layer is repulsive while between different layers 
it is attractive and reduced by a factor s/2X <C 1, where s is the layer spacing and A = X a b 
is the penetration depth for the currents in the layers. As a consequence, the interaction 
of two straight stacks of pancakes at distance p ^> A decays as exp(— p/A) since at long 
distances the attraction and repulsion between the pancakes from both stacks compensates 
almost exactly. Thus, the interaction of two straight stacks is just the usual short-range 
repulsion of Abrikosov vortices. 

The repulsive and attractive interaction of point vortices has a further important conse- 
quence, which to our knowledge has not been pointed out previously, namely, the interaction 
of a distorted pancake stack with a surface which is parallel to the stack. Within the lin- 
ear London theory the condition of zero perpendicular current through a planar specimen 
surface may be satisfied by adding the magnetic fields and currents of image vortex lines 
Each vortex then is attracted to its image since the images have opposite orienta- 
tion (antivortices). However, this short-range attraction applies only when the vortex line 
is perfectly straight. As soon as the vortex line is distorted, the compensation of repulsive 
and attractive terms in the vortex-vortex interaction is no longer ideal. As a consequence, 
randomly distorted vortex lines feel a long-range attraction to the surface. 

In this paper we derive this novel long-range attraction and discuss some of its conse- 
quences. We show that the result of the isolated-layer model exactly coincides with the 
result of anisotropic London theory in the limit A c ^> A, where A c is the penetration depth 
for currents perpendicular to the layers. The long-range attraction is, therefore, also present 
in superconductors with finite anisotropy, where it is proportional to x~ 2 exp(— 2x/A c ) for 
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long distances x between vortex and surface. Since in high-T c superconductors A c typically 
is a macroscopic length, this surface attraction is really of long range. 

We first give a simple physical interpretation of the long-range fluctuation- induced at- 
traction. Assume that only one of the pancakes of a straight stack is displaced by a small 
distance u away from the surface. This local distortion is formally described by adding 
a pancake at the new position x + u and an antipancake at the equilibrium position x, 
which annihilates the original pancake. The same procedure has to be done with the image 
stack situated at the position — x if the surface is at x = (see Fig. [I]). The two pancake- 
antipancake pairs are dipoles with a strength proportional to the displacement u and a 
dipole-dipole interaction energy proportional to u 2 /x 2 , where 1/x 2 is the second derivative 
of the pancake-pancake potential (lnx). Therefore, the distorted vortex line is attracted to 
its image, and thus to the planar surface, by a long-range potential proportional to u 2 /x 2 . 
This long-range interaction is in addition to the short-range interaction of a perfectly straight 
vortex with its image, which is proportional to K (2x/X) oc exp(— 2a;/A), where K (x) is a 
modified Bessel function. 

The interaction of a distorted vortex line with its image may be calculated from the 
interaction of pancakes separated by r mn = (x m - x n ,y m - y n ,z m - z n ) = (x mn , y mn , z mn ) 

EH. 

( 21n(p TOm /£), n = m 

£ ~ C { S / \z mn \\ /p nm \ , (1) 

l"2A exp v"^r) ln l^r)' n ^ m 

with eo = s$Q/(47r/ioA 2 ), p 2 mn = x 2 mn + and z m = ms where s is the layer spacing and 
£ the core radius of the pancake. For z m = z n Eq. (1) applies to all distances p mn > £, but 
for z m 7^ z n p mn ^> A was assumed. The total energy of a vortex line is composed of its 
self-energy and the interaction with its image line of opposite orientation, namely, 

^ m n 

£{x m ~\- X n) y m yni Z>m %n) ] > (2) 

where the sums are over the pancakes of the real vortex in the superconducting half space 



x > 0, thus x m ,x n > 0. For a distorted vortex line parallel to the surface x = we define 
pancake displacements u m = u m (z m ) = (u xm , u ym ) by writing x m = x + u xm and y m = u ym . 
Keeping only the terms quadratic in the u m we obtain from (2) the linear elastic energy of 
the vortex. 

We first consider random and isotropic displacements with ensemble averages (u xm ) = 
(u ym ) = 0, (u xm u xn ) = (uy m Uy n ) = f(\m-n\), (u xm u yn ) = 0. The long-range interaction 
with the image line for 2x ^> A then becomes from (1) and (2) 

where L is the vortex length. This fluctuation- induced interaction is attractive and depends 
on the relative displacements ((u^ — u ) 2 ) over a vortex length of order A. Like a dipole- 
dipole interaction it decreases as 1/x 2 . Since s < A, one may approximate the sum in (3) 
by 



32irfx \ 3 x 2 
where g(z) is the correlation function 



f°° / z \ 

J o dz exp(- -J g(z), (4) 



g(z) = (nz)-u(0)f). (5) 

The integral (4) shows that the long-range interaction does not depend on the layer separa- 
tion s. 

For a more comprehensive calculation we use the general expression for the two-pancake 
interaction |8]||, 

n $^2 d 3 k 1 k 2 

Ho J 8vr 3 1 + k 2 X 2 q 2 1 ' 

where k = (k x , k y , k z ), q = (k x , k y ), p = k z , thus q 2 = k 2 + k 2 and k 2 = q 2 +p 2 , and as above 
r m = (x mn , y mn , z mn ). Since (6) is valid for all r mn , we may obtain the total elastic energy 
of the distorted vortex line E'tot = E se \{ + E int from (2) and (6). Expanding this to quadratic 
terms in u, introducing the Fourier transform 



ufo) = u(p) (7) 

and using s^exp(ip^) = 27r5(p) (for |p| < 7r/s) and f(dp/2ir) |u(p)| 2 = Jdzu(z) 2 = (u 2 )L 
we obtain 

<T)2 r rj^h 

E tot = ~± —\u( P )\ 2 (fs C U + frnt), (8) 
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(9) 



/int = ( l+^A 2 ^ ~ ^ /self ) e2 ^ xX " (10) 

Remarkably, exactly the same result (8)-(10) is obtained from the anisotropic London 
theory in the limit A c —>■ oo. Namely, the interaction energy of two London vortices at 
positions Ti(z) and r 2 (z) is 

Jdri a Jdr 2/ 3 V^ri - r 2 ) , (11) 

where the anisotropic interaction V a p(Ti — r 2 ) (a, (3 = x,y, z) is given in Ref. |TD[ as a Fourier 
integral over V^g(k). In the limit A c — ■> oo, the general expression V^g(k) = ($q//x )(1 + 
/c 2 A 2 )~ 1 5f Q/ g(k) simplifies to a diagonal matrix with g xx = kl/q 2 , g yy = k 2 /q 2 , and g zz = 1 
if z is along the c-axis of the uniaxial superconductor. Inserting this V^g(k) into (11) and 
integrating over the vortex and its image, we reproduce the result (8)-(10) of the pancake 
approach. 

From Eq. (10) one reproduces the long-distance interaction (3), (4). Moreover, from the 
self-interaction (9) the line tension P of an isolated flux line p] or stack of pancakes is 
obtained, 

v° (1 :f 2) - < i2 > 

87r no x z p 2 y 

which determines the linear elastic self-energy of a distorted vortex line, 

EseU=lj^p 2 P(p)\u(p)\ 2 . (13) 



In real space this energy looks similar to Eq. (3), 
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^exp(--|/|J , (14) 
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and since s <C A it can be approximated as 

The similarity of (4) and (15) leads to the following useful relationship. If the correlation 
function (5) increases algebraically, g(z) = const • {zp, one has (still for 2x ^> A), 

£i„t = -(7A74:r 2 )£scif. (16) 

For example, a flux line diffusing in a random pinning potential may exhibit 7 = 1, thus 
E int = — (A 2 /Ax 2 ) .E^if is not a very small correction to E sei{ . 

To obtain the thermal fluctuations one has to consider the pronounced dispersion ( "non- 
locality") of the line tension P oc 1/p 2 (12) at large p ^> A -1 . This means that a single 
displaced vortex "feels" a parabolic potential of curvature p 2 P(p) m $oln(aA/s)/(47r/i A 4 ) 
(with a = 1.16 for s < A) originating from the global interaction with all pancakes within 
a distance of order of A along the stack. Therefore, the self-energy (3) caused by short- 
wavelength fluctuations is simply 

= (U 2) . (17) 

87r/i A 4 w v ; 

Each pancake only weakly feels the usual "local" (non-dispersive) line tension that originates 
from the nearest-neighbor interaction, which here is strongly reduced due to the absence of 
Josephson coupling between the layers. For long tilt wavelengths 2n/q ^> A, however, the 
line tension (12) of the pancake stack is local, P = $Q/(87r/i A 2 ), as can be seen also by 
inserting a uniform tilt — u oc I in (14). 

In terms of the Fourier transform (7) the correlation function g(z) in (4) and (15) may 
be expressed as 

^) = Z/2| (|G(P)|2)[1 " C0S(PZ)] - (18) 
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Writing the self energy (11) for an isolated vortex of finite length L as a sum over discrete 
modes with wave vectors p n and amplitudes u n and ascribing to each mode the average 
thermal energy ksT (since u ra has two components), one obtains (|u(p)| 2 ) and with (18), 

327r/i A 4 [dp l-cos(pz) 

9{z) = kBT hi+p^y (19) 

The function g(z) is practically constant in the interval relevant in Eqs. (4) and (15). 
One has g(0) = and for \z\ ^> s one finds g(z) ~ c\ + c-2\z\j\ with ci <C c\. For the 
example \/s = 100 some values of the integral in (19), times A, are 0, 4.71, 5.07, 5.16, 5.25, 
5.29, 5.42, 5.53, 5.68, 5.84, 6.10, 6.85 for z/s = 0, 1, 2, 3, 4, 5, 10, 20, 50, 100, 200, 500. The 
constancy of g at z up to a few A means that the thermal fluctuations of the pancakes of one 
vortex line are nearly uncorrelated at small wavelengths P,|12||, and thus the cutoff length of 



several s, or the effective cutoff \q\ < a/s, required to obtain expressions like Eq. (17) is not 
crucial. This "elastic independence" of the pancakes originates from the strong dispersion 
of the line tension P(p) (12). For \z\ ^> 27rA and for a usual string with non-dispersive P 
the thermal random walk yields g(z) oc \z\. 

The expression for g(z) (19) was derived from the line energy E 8e if (13), which does not 
depend on the vortex distance x from the surface. The correct calculation, however, has to 
consider the total energy E se if + E int . Since the interaction with the surface E int depends on 
x, and the equipartion theorem requires constant energy ksT per mode, we find that the 
fluctuation amplitude depends on x. Explicitly we obtain from (4) and (15) for a vortex line 
at x ^> A/2 the short-wavelength thermal fluctuations, i.e., the value g(z) at s < \z\ < A, 

2 _ 2 A 2 k B T r ^ 1 I 
(U } ~ \n(a\/s) "17 y + 4^2 \ n (a\/s) i ' { } 

So we have the interesting result that the correction to the thermal fluctuations decreases 
away from the surface only as a power law. This increase of (u 2 ), which in the total energy 
exactly compensates the spatial dependence of the term E int , originates from the softening 
of the flux-line lattice near the surface. 

Considering u x and u y separately, we find that the surface-caused softening of a vortex 
line is indeed isotropic at large x but becomes anisotropic for x < A. The general expression 



obtained from (2) and (6) shows that to a good approximation the contributions of each 
displaced pancake to the elastic energy add independently. Each pancake contribution is 
composed of the interaction between six objects, namely, the displaced pancake and the hole 
(antipancake) it leaves in the vortex line, their images in the same layer, the undisplaced 
straight vortex line, and its image. Explicitely one has 

E 



,0 ' " 87T/i A 4 



a(x)(u x ) + b(x)(u y ) 
aX X 2 rr „/2x 



a(x)=ln^-^-^m 

l2 



Ax 2 * V A 1 



, . . , aX A A r . /2x\ , , 

with a = 1.16. For large x ^> A/2 one has a = b = ln(aA/s) — A 2 /4x 2 , and for x <C A/2, 
a = b — X 2 /2x 2 , b = ln(aA/s), since the derivatives are K' (r) = — l/r and K{j(r) = 1/r 2 for 
r <C 1. The thermal fluctuations 

2 2k B TX 2 2 2k B TX 2 

K) = 1^T' {Uy) = ^bW (22) 

become thus asymetric near the surface, with (u 2 )/ (u 2 ) = b/a > 1. Note that the curvature 
a in (21) turns negative at small x. To check the stability of a pancake stack close to 
the surface and to obtain the correct pictures of penetration and exit, one has to add the 
(always positive) stray- field energy [0,0 and the restoring force of the current density j(x) 
that originates from the Meissner surface current and from other vortices, replacing a by 
a + (47r/i A 4 / <: I ) o)j '{ x )- Work in this direction is under way. 

Large fluctuations and local tilt of flux lines can be caused by pinning. Random pinning 
forces on a single flux line would cause square displacements diverging proportional to the 
number of forces. More realistic random pinning potentials lead to finite vortex displacements 
P1J15H- For a crude estimate assume that random pins of density n p are so strong that the 
vortex line wanders an average distance squared of order (4n p s) _1 as it passes to the next 
layer. This yields g(z) w z/(n p s 2 ), and the interaction (4) with the surface becomes 

E ® lL ^ E (2^ 

int " 32n n X*n p s 2 x 2 ~ 4x 2 self ' { 6) 
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The resulting long-range force —dE- mt /dx in principle may drive the vortex line to the surface. 
It has to be added to geometric forces which are exerted on the vortex line by surface 
screening currents when the specimen has a non-ellipsoidal cross section and is not infinitely 



thick 

The elastic energy (8)-(10) of a single flux line may be generalized to the elastic energy 
of the distorted flux-line lattice in the half space x > by replacing in it |u(<?)| 2 by the 
three-dimensional Fourier transform |u(k)| 2 of the lattice displacements ||16|| . The complete 



expression for E to t slightly differs from our approximation (8), which applies if (it 2 .) = 
(Uy) = |(u 2 ) and (u x u y ) = 0. The presence of a flux-line lattice of density B/Q increases 
the fluctuation-induced attraction to the surface (3) when the displacements in the same 
layer are correlated over a radius p corr > (Qq/B) 1 / 2 . The interaction per vortex will then 
increase by a factor of order of 7rp 2 orr .B/$o- 

Similar expressions for the elastic energy of isotropic and anisotropic vortex lattices 



near a planar surface are derived in Refs. [13||T4[] . In isotropic superconductors a long-range 



attraction of the type (3) is absent, and in anisotropic superconductors the attraction is 
reduced by a factor of exp(— x/X c ). Formally, the long range is due to the non-cutoff factor 
1/q 2 in Eq. (10). 

In conclusion, we have shown that in extremely anisotropic layered superconductors there 
is a long-range interaction between the sample surface and a distorted vortex line parallel 
to the surface and to the c-axis. This interaction causes a spatial variation of the thermal 
fluctuations even at distances much larger than the in-plane London penetration depth, 



which might affect the melting process of the vortex lattice [19] and the evaporation of 



vortex lines into independent pancake vortices []3,|l2 . In the presence of sufficiently strong 



random pinning, a fluctuation-induced long-range force attracts the distorted vortex line to 
the surface. This additional force means a bulk current density far inside the superconductor 
and may lead to flux creep towards the surface. 

We acknowledge support from the German-Israeli Foundation for Research and Devel- 
opment, Grant # 1-300-101.07/93. 
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FIGURES 

FIG. 1. Left: A distorted vortex line and its image composed of pancakes (|) and antipancakes 
{[). Right: The two dipoles generated by the displacement cause a long-range attraction between 
the distorted vortex and the surface (indicated by a vertical bold line). 
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